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Abstract 



We show that the radial Teukolsky equation (in the frequency domain) with 
sources that extend to infinity has well-behaved solutions. To prove that, we 
follow Poisson approach to regularize the non-rotating hole, and extend it to 
the rotating case. To do so we use the Chandrasekhar transformation among 
the Teukolsky and Regge- Wheeler-like equations, and express the integrals 
over the source in terms of solutions to the homogeneous Regge- Wheeler-like 
equation, to finally regularize the resulting integral. We then discuss the 
applicability of these results. 
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I. INTRODUCTION 



With the approach to the beginning of a new era in gravitational wave detectors, a 
renewed interest in the theoretical study of astrophysical sources of gravitational radia- 
tion arouse. Several groups of research are now reviewing computations carried out twenty 
or more years ago. One of the most remarkable "old" results, is the description of the 
gravitational perturbations about black holes. In particular, for the most plausible astro- 
physical scenario, i.e. the Kerr (rotating) hole, perturbations are compactly described by 
the Teukolsky equation Among the new theoretical results, stands out the close limit 
approximation [Q, which approximates the collision of two black holes by a single perturbed 
one. The subsequent extensions to second order perturbations [^, and to moving (head-on) 
holes IQ, confirmed the success of the approximation. Also recently the collision of a 
small and a big black hole by perturbative methods have been revisited in order to incorpo- 
rate non-vanishing initial data. In both approaches only perturbations around Schwarzschild 
(non-rotating) black holes have been studied by use of the Zerilli-Moncrief equation: A to- 
tally well-defined description. The success of these techniques encourages now the extention 
of them to a more realistic rotating background. In the way to obtain such generalization 
lie some analytic and conceptual problems. One of them, the regularization of the Teukol- 
sky equation, have been recently clarified by Poisson ||^. Apparently, it was first stressed 
by Detweiler and Szedenits 0, that when one tries to find a solution (by straightforward 
application of the Green method) to the radial Teukolsky equation with a source term that 
extends to infinity, one obtains a divergent result. Not much later, Sasaki and Nakamura 
IP realized that applying a Chandrasekhar transformation to the radial non- homogeneous 
(including the source term) Teukolsky equation one could end up with a Regge- Wheeler-like 
equation (with a potential term that reduces to the Regge- Wheeler potential in the a — > 
limit), that had a real, short range potential (as opposed to the original Teukolsky one), and 
generated sensible finite results for extending-to-infinity sources. In this work we generalize 
Poisson regularization procedure (that he applied to the non-rotating version of the Teukol- 
sky equation [l^), to the rotating case, which is the relevant case for further applications, 
since otherwise one would preferably use the Zerilli-Moncrief formalism. 

The paper is organized as follows: In Sec. II we review the Teukolsky equation in its 
time and frequency domains as well as the asymptotic behavior of its homogeneous solutions. 
Section III is devoted to recall the main results of the Chandrasekhar transformation that 
maps solutions to the radial Teukolsky equation into solutions to a Regge- Wheeler-like 
equation. In Sec. IV we formally follow Poisson regularization procedure for the a 7^ 
case and in Sec. V we explicitly apply it to the problem of a particle infalling along the 
symmetry axis of a Kerr hole {6 = 0). Finally, a brief discussion of the applicability of these 
results is given. 

Notation and conventions: The sign conventions for the metric and curvature tensors are 
those of Misner, Thorne and Wheeler That is, the metric signature is (— + + +); 
the Riemann tensor is defined by the equation -R"^^^ = T"^^ 5 - r"/?^,^ + P'^^^P'^^ - P'^^^P"^^ 
and the Ricci tensor by equation R^^, = R^^^u- The tortoise coordinate r* is defined so 
that d/dr* = f{r)d/dr, where /(r) is specified in Sec. II. A. A prime indicates differentiation 
with respect to the radial Boyer-Lindquist coordinate r. An overbar represents complex 
conjugation. 
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II. THE TEUKOLSKY EQUATION 



In 1972 Teukolsky |]T2|, using the Newman-Penrose formalism, succeeded in desentan- 
gling the perturbations of the Kerr metric, and wrote a master equation for fields of arbitrary 
spin (s = for the scalar field, s = ±1 for the electromagnetic field, s = ±2 for the gravita- 
tional field, etc). In Boyer-Lindquist coordinates {t, r, 9, ip) this equation, for pure "outgoing" 
(s = —2) gravitational perturbations, takes the form 
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- (4 cot' ^ + 2) = 47rST, 



where M is the mass of the black hole, a its angular momentum per unit mass, S 



r' + a' cos' 9, and A 



2Mr + a\ 



The field \1/ = \l/()f:, r, ^, Lp) is the outgoing radiative part of the Weyl tensor 



(2.2) 



i. e. one of the two gauge invariant components of the perturbed Weyl tensor along the 
Kinnersley tetrad vectors n° and m", with p being a spin coefficient. The explicit expression 
of the tetrad vectors and of the spin coefficients are given in Table |, and the form of the 
source term, T = T(t, r, 6*, y?), in equation ( |2.1| ) is given in Section V, Eq. (p.l|). 

Equation ( p.lj) can be separated by mode decomposition of field and source term in the 
following way 
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The radial functions RirnA^{r) satisfy the following equation ]TB[ ] 

o d ( 1 d 



dr \Adr) ''^^""^ 
with a complex effective potential of the form 



f/(r) 



A 



[cjV + 2zu;o-A' - {Aiua + A) A], 



(2.5) 



(2.6) 



where a = (r' + a') — am/uj, X = E — 2 + ^^'(a' — am/uj). E is a, separation constant [|17| , 
that in the case a = is given by /(/ + 1) and the source term is given in Section V, Eq. 
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The angular functions Sf_2dmi^) called spin- weighted spheroidal harmonics and sat- 
isfy an eigenvalue equation of the form 

J 1 ^ ( ■ fm'^ + 4- 4:171 COS 9 \ 

+ a'u;'cos'e + Aaujcose'^S^,^i^je) = -ES\i^^{e). (2.7) 

In the case atu = they reduce to the spin- weighted spherical harmonics, i.e. imi^^^^ = 
0) = .2Yi^=oiO). 



A. The boundary conditions for the Teukolsky equation 

The inhomogeneous radial equation can be solved with the physical requirement that 
gravitational waves must be purely ingoing at the black hole horizon, located at r_|_ = M + 
y/M"^ — a?, and be purely outgoing at infinity. In mathematical language these statements 
translate into 

r.H . *^ _ / Pe-'^'* r^r+ (r* ^ -oo) . . 

H,rau^\ ) \ Q^{iojr)-^e-'^''' + Q^"* {iurfe'^""' r ^ oo {r* ^ oo) ' ^ ^ 



and 



l«u;r)^e*'^'^* r ^ cxd ir" ^ ool 



Rl,7nu;{r ) ~ I pin p^-iur* _^ pout ^iur* j. j.^ _^ ' (2-9) 

where Q!f and , and P^*^* are constants (independent of r) and /(r) = A/(r^ + — 
am/uj). The coordinate r*, as defined by Chandrasekhar is 

dr /A\ M2-am/(2tu), (r-r,\ 

-r + Mln^ + ,,,, ^\ M n i , (2.10) 



where r± = M ± V — a^. Note that this is not the standard tortoise coordinate, which is 
instead defined by / (ir(r^ + a^)/A. The reason of our choice for this non-standard definition 
of r* (defined in Eq. ( p.lOj )) lies in the simplification of the forthcoming equations. From 
the numerical point of view this choice should not be especially inconvenient. 
The wronskian of these solutions is given by 

WTiR"", R°°) = (RhR'oo - R^R'h) = -2iuj''Ql^A. (2.11) 



of I^T, i.e. -22^3 gm_ 



Note that the "conserved Wronskian" as considered in Refs. ||7|,pr3| is only the constant factor 
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III. CHANDRASEKHAR TRANSFORMATIONS 



It is called a Chandrasekhar transformation |I^ the first order differential operator that 
carries a solution to a Regge- Wheeler like equation (with a real potential) into a solution to 
the homogeneous Teukolsky equation. More explicitly 

R^'^{r) = x^'°°CX^'°°(r), (3.1) 

where x^'°° are constants of normalization (see the Appendix) and C = Gconv{r) + Gdiv{r)C 
with C = d/dr* + iu is the first order operator. The form of the functions Qconvi'f') and 
Gdivi'f') (made explicit in the following subsection) are determined by the requirement that 
the potential term, V{r) of the resulting Regger- Wheeler- like equation 

CCX = V{r)X, (3.2) 

be real. 

The asymptotic behavior of the solutions to this equation can be described in terms of 
two solutions representing purely ingoing gravitational radiation at the black-hole horizon, 
denoted by X^{r) and gravitational waves purely outgoing at infinity, denoted by X°°{r). 
Its asymptotic behavior can be characterized by 

_ / e-^'^'-* r ^ r+ (r* ^ -oo) , . 



Qrw e-"^" + Qrw e*^' r ^ oo (r* ^ oo) 



and 



[iRW^ + I RW^ r ^ r+ [r ^ — ooj 

where Q'r^ and Qrw^ '^rw and are constants (independent of r), which can be related 
to the corresponding Teukolsky constants, Q™, QI^"*, and P^*^*. 
The wronskian of these solutions is given by 

^RW - - X-^ = 2^u^Ql^ . (3.5) 



A. Explicit form 



In 1975 Chandrasekhar was able to reduce the Teukolsky's radial equation governing 
the general, non-axisymmetric, ingoing (s = -|-2) gravitational perturbation of the Kerr black 
hole to a one- dimensional wave equation of the form (p.2|). For the sake of future numerical 



convenience, we re-express, here, his results for the case of outgoing radiation (s = —2). 



The functions that characterize the transformation (3.1) are the following 
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where /(r) = A/(r^ + — am/uj). q and p are functions of r, and h and are constants, 
all as defined below 



q{r) = {\- 2)or2 + 3or(r2 - a^) - Sr^A , 
p{r) = q'A — qA' , 

b = ±3 (a'^ — am/u?j , 

A; = ± {36M^ - 2(A - 2) [(a^ - am/cj)2(5A - 4) - 12a^] + 26A(A - 2)} 
The potential in the Regge- Wheeler equation is then given by 



(3.7) 



-1/2 



V{r) = A 



A A(A 

-b— + ^ 



2) {p/a - kA) {kq - bp/a) 



q + bA (q + bA) {q-bAf 



(3.8) 



Note that since b and k can be chosen with either sign, Eq. ( p.8|) yields, in general, four 
different potentials. For the choice of the positive sign in the expression for k, one obtains 
two Zerilli-like potentials (depending on the sign of b), which merge to the known Zerilli 
potential (for polar or even parity perturbations) in the limit of a ^ 0. For the choice of 
the negative k, one obtains two possible Regge- Wheeler-like potentials, which reduces to the 
known Regge- Wheeler potential (for axial or odd parity perturbation) in the limit of a — * 0. 
Although different, the four potential have the same asymptotic features of vanishing at the 
horizon and of falling off as (9(r~^) at infinity. 



IV. A GENERAL PROCEDURE OF REGULARIZATION 



To regularize the solution of the radial Teukolsky equation with sources that extend to 
infinity, we can proceed along the same lines of Poisson |^ and generalize his procedure for 
the rotating case. It is also convenient to use here a notation parallel to the one used for 
the regularization of the Bardeen-Press equation. For instance, we can express the radial 
source term (see Eq. (|5.6| )) in the following form 



Tjr) 
A2 



Gg{r)e 



iLot{r) 



(4.1) 



An explicit form for T(r) or equivalently for the constant G and the function g{r), is given in 
section V. For the purpose of this section its sufficient to anticipate that A~^T(r) ~ (9(r"^/^) 
as r — > oo. 

The general solution to the inhomogeneous Teukolsky equation can be expressed as 

GAj 



R(r) = ^|i?°°(r) A + g{f)R^{f)e''^'^'Uf 



B 



^^(f)i?°°(f)e*^*(") df 



(4.2) 



where A, 5, c, and d are constants to be determined by the boundary conditions. Note that 
if we want to impose here the boundary conditions (|2.8|) and (|2.9|) , i.e. take A = B = 



6 



c = r+, and d = oo, the integrals do not converge when the superior hmit goes to infinity. 
We thus have first to regularize them for only then to impose the boundary conditions. Such 
process of regularization leads to Eq. ( [4.12| ) below. 

We apply the Chandrasekhar transformation C given in Eq. ( p.l|) to the following integrals 



r_ff,oo_ 



G / g{r)R 



H,ooi 



re 



where 



I g{r) (g^conv + Qa,X) X^'°°e*-*(^)rfr 

X G(^I conv ~l~ -^div) ; 

. = G'|(7(r)6;,onve^-*MX^'°°(r)dr, 



(4.3) 



(4.4) 



is a convergent integral since the integrand goes as 0{r ^/^) for r ^ oo (from Eq. ( ^.6p , 
Gconv ^ 0{r)). In this way we can restrict our attention to the divergent integral 



div 



G J (?(r)^?dive'"*(")/:X^'°°(r)dr 



(4.5) 



where now the integrand is 0{r ' ) as r — > oo (from Eq. ( p.6| ), Qdw ~ 0{r )). 

As noted by Poisson 0, Jjiv can be regularized by introducing a function h{r) and 
integrating by parts i. e. writing it in the following form 



dr - he'^^'CX 



boundaries 



/div - he"^'CX 



(4.6) 



boundaries 

where h{r) is a function to be determined. Using the following property for the operator C 

1 



dr 



fir] 



{V{r)+iujC)X 



the new divergent integral becomes 
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THi 



div 



dh . { dt 1 ■ , , ^ / X ^ 
— + tuj( — + -]h + Gg (r ) ^div 
dr \dr j, 



dr 



0) 



(4.7) 



(4.8) 



(4.9) 



Here dt/dr is the equation of the trajectory, explicitly given in Section V, and the function 
h{r) must be chosen so that /div converges when r oo. The natural and elegant choice of 
Tdiv = removes completely the divergence in Jdiv as in the Schwarzschild case studied by 
Poisson f^, 
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h{r, a = 0) = 8M' 



2 1 + X + 2iMujx^ 
1 + x 



X = {2M/r)-'/\ 



(4.10) 



but the power counting shows that it is enough to remove all the dependence up to the 
first power in a, i.e. terms like (a/r)"^ or higher produce convergent integrals. Although is 
possible to solve formally the consequent first order differential equation for h{r), it's a very 
difficult task to find the explicit form of the particular solution h{r), with a 7^ 0. To handle 
this technical difficulty, one can, for instance, numerically integrate this differential equation 
along with the homogeneous Regge- Wheeler-like equation to obtain the final solution ( [4.21) . 

Using now Eqs. (^]3|), ( [4.4| ), ([4.6|), and ( [4.8| ) we can express Eqs. ( [4.3| ) in the following 
regularized form 



rH,ao 



X 



H,oo 



fir) 



- X 



H,oo 



boundaries 



(4.11) 



In summary, after the imposition of the boundary conditions ( p.8| )-( pl9D , the general 
solution ([4.2|) to the radial Teukolsky equation can be written in the following regularized 
form 



X^i?~(r) 




[ (G^7(r)^conv + jrh'ir) ) X^(r~)e^-*(^-) df 



fir 



h{r 



/(r 



(4.12) 



+ x"x^ 



where the boundary term has been simplified by expressing R^''^{r) in terms of X^'°°(r) 
by the transformation ( p.l[ ), and using the following commutation relation [C,X^'°°] = 
f{r)dX"'°°/dr. 

Note that the integrals are now not only convergent, but also expressed in terms of the 
X{r), solutions to the Regge- Wheeler-like equation, i.e. a one- dimensional equation with a 
real and short range potential. 

To have an expression for the regularized integrals we have to choose a Chandrasekhar 
transformation that gives an explicit Qconv and Qdiv, and produces a potential V{r). The 
completion of the solution is obtained by determining the normalization constants fhat 
relate the behavior of the R^'°° and the X^'°°. In the Appendix we describe how to find 
them. 



V. THE SOURCE TERM 

In order to see explicitly how the regularization works, we must construct the source term 
of the Teukolsky equation for the specific case of a particle of rest mass rriQ released from 
rest at infinity and falling with zero angular momentum along the polar axis of a Kerr black 
hole of mass M and angular momentum J = aM. Note that this choice is only dictated 
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by the simplicity of the resulting source term, but the procedure, as explained in Sec. IV, 
works for any source extending to infinity. 

The source term in the Teukolsky equation T = T(t,r,9,ip), can be constructed 

from the following expression [|12| 



T = 2p-^| (-A + 37-7 + 4/1 



S-2f + 2a] Tnm - (-A + 27 - 27 + /i) 



+ -5 - f + 5/3 + 7a 



-A + 27 + 2/i T, 



-5 -f + 2(3 + 2a]T, 



(5.1) 



where the spin coefficients, operators (denoted by an overhat) and Kinnersley tetrad are 
explicitly given in Table |. 

To begin with, we replace into it the following non-vanishing NP components of the 
stress-energy momentum tensor 



-'-mm — ■'-fj.u'"' '"' 



-'-nm — -'-fj.u"' '"' 



a^p^ sin^ 9 



iap sin 6* 
2v^S 



Tu(r^ + a^)-TrtA 



(5.2) 



T —T n^n" 

nn — fiv''' '<' 



4S2 



Tttir' + ay + TrrA' - 2Trt{r' + a')A 



We then consider the stress-energy-momentum tensor for the trajectory specified above 



T^" = ^Kr)5{r - R{t)) 5 {cos 9 - 1)5(<^), 



(5.3) 



where /(r) = A/(r^ + a^). 

The only non vanishing components of the four-velocity are 

CLT 

u'^ ^ = -(1 - hr))-"\ 



The equation of the trajectory is 

dt _ 
dr 



-f{r)-\l-~f{r))-'l\ 



(5.4) 



(5.5) 



Finally, the source term of the radial equation (|2.5|) can be computed from the inverse 
Fourier transform of Eq. (|2.4| ) , 



2 J dtd'fd{cos9)e'^'''-""^^J:TS^2,i,m{0)- 



(5.6) 



Using symmetry arguments (m = for the particular trajectory we have considered here) 
we then obtain (in the notation of Eq. ( [4.1|) ) 



9[r) 



[r — la) 
4r2A 



[(2-/)(l-/)-V2_2 



(5.7) 
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and 

G = -16moN%i.Q{auj) , (5.8) 

where the constant N^i.Q{auj) is defined by the hmiting behavior of spin weighted the 
spheroidal harmonics |T9| 

lim5-(^)=^l-+^liV-o(au;). (5.9) 

In the case au = one has 



V An 



VI. SUMMARY AND CONCLUSION 



The procedure for regularizing the Teukolsky equation with source terms that extend to 
infinity (infalling particle) in the frequency domain was straightforward but lengthy. The 
finite result is given in Eq. ( [4.12| ), where G, Wt, and Wrw are given in Eqs. (|578|), ( p. Ill) 



and i WM respectively, and in the Appendix. R^^'^ir) and X^'°°(r) are solutions to 

the radial Teukolsky, Eq. (|2.5|) and Regge- Wheeler, Eq. (p.2|), homogeneous equations. Qc 



7conv 



and Qcorw are given in (p.6|) for the specific transformation chosen by Chandrasekhar [|I4 
g{r) is given in Eq. ( |5.7| ), and bring all the dependence on the source, and hence on the 
particular trajectory of the particle one wants to study. Finally, the regularizing function 
h{r) can be determined by solving the first order differential equation ( |4.9| ). We repeat here 
that although we have studied explicitly only the case of the particle infalling from rest at 
infinity, it is obvious that the regularization method applies to any type of trajectory (the 
final expressions for will be more complicated though). 

The method of regularization used in this paper has been chosen to parallel that of 
Poisson [^. It is clear that one could have used a different procedure to regularize directly 
the Green function, as one does in field theory. The final result, of course, has to be 
independent of the method chosen. It remains to be understood what is explicitly the Green 
function of the radial Teukolsky equation in the case of sources that extend to infinity. We 
would like to come back on this issue in the near future. 

This long process to determine the regularized solution, although valid in general, does 
not seem very practical at the moment of finding a numerical solution to a given astrophysical 
problem. In this sense, if one chooses to work in the frequency domain, the approach by 
Nakamura and Sasaki ^ , that find a Chandrasekhar transformation for both, the waveform 
and a generic source that produce finite results, appears more appropriate. In situations 
where one has to work in the time domain, like when initial data are not negligible, one 
has to use the Teukolsky equation. It is not yet known if in this case one has to perform 
again a regularization of the results. Besides, one can actually implement a procedure in 
order to obtain numerical results and compare them with those obtained by solving the 
Sasaki-Nakamura equation (for the particle infalling along the symmetry axis). This will be 
the subject of research in our next paper. 
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APPENDIX A: NORMALIZATION CONSTANTS 

To find the two normahzation constants ^^id we need to look at the asymptotic 
behavior of the solutions to the wave equation (p.2|), given in Eqs. (p.3|)-(p^. The asymp- 
totic behavior at infinity do not pose us any problem and the result of = 1/4, found by 
Poisson is valid here as well as in the a 7^ case. This is so because the effects of rotation 
on R°°{r) vanish with a higher negative power of r than the corresponding nonrotational 
ones. 

To study the asymptotic behavior of the solutions near the horizon located at r = r+, let 
X{r) = y(/)e~*'^''*. A simple change of variables tell us that Y{f) satisfies a second order 
differential equation 

ff'Yff + [ff" + //'(/' - 2^^)] Yf - yU)Y = 0. (Al) 

The solution to this differential equation can be given as a series of the form Y{f) = 
1 + o-nf"'- The first two coefficients are then found to be 

^(^^ 

CLi 



a2 



/('o)(/('o)-2za;)' 

^{2)-ai [/{0)+2/('i)(/('o)-M-l^(l) 



4/(o)(/(o)-^^) 



(A2) 



where we have used the notation V{f)/ f = V(i) + V(2)/ + ... (V(o) = since V{r), given in 
Eq. (|3.8|), vanishes at r = r+). And the derivatives of / with respect to r are 



/(o) 



2VM2 - a2 



f'm 




1 4rj 



where the subindex "+" means quantities evaluated at r+. 

We can now apply the Chandrasekhar transformation to and determine the normal- 
ization constant such that = x^CX^ reproduces the asymptotic behavior (|2.8|) . We 
then have asymptotically close to the horizon 
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= X^'e-^-'-*! (6;,„(o) + g^MDf + gd^vi2)f + 0{f)) {aj + 2a2p + 0{f)) f 

+ (6^co„.(0) + e^con.d)/ + Gconvi2)f + 0{f)) (l + + a^f + |. (A4) 

Comparison order by order gives 

conv{l) 

, (A5) 

on the horizon, and 

= \^Gconv{2) + fll Gconv{l) + Gdiv{l)f(o) + Gdiv{0)f{l) + 2a2^(iii;(0) /(O) } • (^6) 

The exphcit computation of the series developments of G and V is tedious but straight- 
forward. Direct use of Eq. ( |3.6| ) gives 



Gconv(l) ~ 


= luja^ V(i) , 


Gdiv{0) = 






-1/2 . / 

lua^ f {q+ 


Gdiv{l) = 


2 I 



^ /(O) ^+ 



r3 T/,,,\ 



(A7) 
(A8) 



^(1) = — 



—^ + -^{2ru;a^-^\)\. (A9) 
^+ /(o) 



(AlO) 



Now, from the expression of the potential ( p.8| ) we can obtain 



^(2) = TfV^(i)--^-A(A-2)-± 



AVg3 

p+{k a+q+ - hp+) f 



+ ^^i^^^^^^^ ( 2^g, + 3g; - 6A; V 



where the subindex "+" means quantities evaluated at r+, and p = q'A — qA'. 
When a = expressions notably simplify and give 

H ^ (l-2^Mu;)(l- 4zM^) 

^ iMuj[{l-l)l{l + l){l + 2)-12iMuj]' ^ ' 
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TABLES 



TABLE L 



Kinnersley tetrad 


operators 


spin coefficients 




= [r2 + a2,-A,0,a]/(2E) 


A = n'^da 


p = 1/ {r — ia cos 6) , r = 


iapp sm6/y/2 


= p/ ^/2[iasui 6, 0, sin 9] 


6 = m°'da 


(3 = pcote/2/V2 , p 


= -p2pA/2 


P = [(r2 + a2)/A,l,0,a/A] 


D = r^a 


a = p{cote + 2iapsme)/2/V2 , 7 


= -p- ppA'/4 
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